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Abstract 

We study the ground state solutions of the Dirac-Fock model in the case of weak 
electronic repulsion, using bifurcation theory. They are solutions of a min-max problem. 
Then we investigate a max-min problem coming from the electron-positron field theory 
of Bach-Barbaroux-Helffer-Siedentop. We show that given a radially symmetric nuclear 
charge, the ground state of Dirac-Fock solves this max-min problem for certain numbers 
of electrons. But we also exhibit a situation in which the max-min level does not 
correspond to a solution of the Dirac-Fock equations together with its associated self- 
consistent projector. 
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1 - Introduction. 

The electrons in heavy atoms experience important relativistic effects. In compu- 
tational chemistry, the Dirac-Fock (DF) model T, or the more accurate multiconfig- 
uration Dirac-Fock model (2j, take these effects into account. These models are built 
on a multi-particle Hamiltonian which is in principle not physically meaningful, and 
whose essential spectrum is the whole real line. But they seem to function very well in 
practice, since approximate bound state solutions arc found and numerical computa- 
tions are done and yield results in quite good agreement with experimental data (see 
e.g. IS]). Rigorous existence results for solutions of the DF equations can be found in 
jll and 15'. An important open question is to find a satisfactory physical justification 
for the DF model. 

It is well known that the correct theory including quantum and relativistic effects is 
quantum electrodynamics (QED). However, this theory leads to divergence problems, 
that are only solved in perturbative situations. But the QED equations in heavy 
atoms are nonperturbative in nature, and attacking them directly seems a formidable 
task. Instead, one can try to derive approximate models from QED, that would be 
adapted to this case. The hope is to show that the Dirac-Fock model, or a refined 



version of it, is one of them. Several attempts have been made in this direction (see 
^ Xi &i .9 and the references therein). Mittleman 0, in particular, derived the DF 
equations with "self-consistent projector" from a variational procedure applied to a 
QED Hamiltonian in Fock space, followed by the standard Hartree-Fock approximation. 
More precisely, let H'^ be the free Dirac Hamiltonian, and a perturbation. We denote 
A+(il) = X(o,oo)(^^^ + The electronic space is the range 7i+(ri) of this projector. 
If one computes the QED energy of Slater determinants of N wave functions in this 
electronic space, one obtains the DF energy functional restricted to (7i^"(ri))^. Let ^si 
be a minimizer of the DF energy in the projected space (7i^(ri))^ under normalization 
constraints. It satisfies the projected DF equations, with projector A"*" (il). Let E{n) := 
f (^'si). Mittleman showed (by formal arguments) that the stationarity of iJ(Sl) with 
respect to ft implies that A+(ri) coincides, on the occupied orbitals, with the self- 
consistent projector associated to the mean-field Hartree-Fock Hamiltonian created by 
From this he infers ([0], page 1171) : "Hence, 51 is the Hartree-Fock potential 
when the Hartree-Fock approximation is made for the wave function". 

Recently rigorous mathematical results have been obtained in a series of papers by 
Bach et al. and Barbaroux et al. |10[ 1111 112| on a Hartree-Fock type model involving 
electrons and positrons. This model (that we will call EP) is related to the works of 
Chaix-Iracane [S] and Chaix-Iracane-Lions Note, however, that in ^1^2^] the 
vacuum polarisation is neglected, contrary to the Chaix-Iracane approach. In |1U| . in 
the case of the vacuum, a max-min procedure in the spirit of Mittelman's work is intro- 
duced. In ^21 J in the case of A^-electron atoms, it is shown that critical pairs (7, P"'') 
of the electron-positron Hartree-Fock energy Sep give solutions of the self-consistent 
DF equations. This result is an important step towards a rigorous justification of Mit- 
tleman's ideas. All this suggests, in the case of A^-electrons atoms, to maximize the 
minimum E{il) with respect to 51. It is natural to expect that this max-min procedure 
gives solutions of the DF equations, the maximizing projector being the positive pro- 
jector of the self-consistent Hartree-Fock Hamiltonian. We call this belief (expressed 
here in rather imprecise terms) "Conjecture M" . 

In [21 and when analyzing the nonrelativistic limit of the DF equations, Es- 
teban and Sere derived various equivalent variational problems having as solution an 
"electronic" ground state for the DF equations. Among them, one can find min-max 
and max-min principles. But these principles are nonlinear, and do not solve Conjecture 
M. 

In this paper we try to give a precise formulation of Conjecture M in the spirit 
of Mittleman's ideas and to see if it holds true or not, in the limit case of small 
interactions between electrons. We prove that in this perturbative regime, given a 
radially symmetric nuclear potential. Conjecture M may hold or not depending on the 
number of electrons. The type of ions which are covered by our study are those in which 
the number of electrons is much smaller than the number of protons in the nucleus, 
with, additionally, c (the speed of light) very large. 

The paper is organized as follows : in §2 we introduce the notations and state our 
main results (Theorems 1^ and II 1(1 . Sections 3 and 4 contain the detailed proofs. 

2 - Notations and main results. 

In the whole paper we choose a system of units in which Planck's constant, h , and 
the mass of the electron are equal to 1 and Ze^ = 47reo, where Z is the number of 
protons in the nucleus. In this system of units, the Dirac Hamiltonian can be written 
as 

H" = -ica-V + c^l3, (1) 
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where c > is the speed of hght '/'—^q '^l^''^^ '^^ ~ ct 'o^ ) 

and the ais are the PauU matrices. The operator H'^ acts on 4-spinors, i.e. functions 
from R3 to C"^, and it is self-adjoint in L'^{M.^,C'^), with domain H'^iR^X^) and form- 
domain ffi/2(K3,C''). Its spectrum is the set (-00, -c^] U [c^, +00). 

In this paper, the charge density of the nucleus will be a smooth, radial and com- 
pactly supported nonnegative function n, with J n = 1, since in our system of units 
Ze^ = AneQ. The corresponding Coulomb potential is V := — n * (l/|a;|). Then 
y : M'^ ^ (— cx), 0) is a smooth negative radially symmetric potential such that 

— -p-r < V(x) < (Vx) , |a;| V(x) ~ —1 for \x\ large enough . 
\x\ 

Note that the smoothness condition on V is only used in step 3 of the proof of Propo- 
sition El Actually we believe that this condition can be removed. 

It is well known that H'^ + V is essentially self-adjoint and for c > 1, the spectrum 
of this operator is as follows: 



(T(iJ" + V)^ (-00, -c^] U {XI, A^, . . . } U 
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with < A5; < A2 < . . . and hm = c . 

Finally define the spectral subspaces = Ker(_ff + V — X^l) and let N^ denote 
tVI^'s dimension. 

Since the potential is radial, it is well known that the eigenvalues A^ are degenerate 
(see e.g. ^1]). For completeness, let us explain this in some detail. To any A G SU{2) 
is associated a unique rotation Ra £ 5*0(3) such that Va; £ M'^, {RAx)-a = A{x-a)A~^ , 
where a = (cti, cr2, fa). This map is a morphism of Lie groups. It is onto, and its kernel 
is {/, — /}. It leads to a natural unitary representation • of SU{2) in the Hilbert spaces 
of 2-spinors L^iS^,C^) and L^{R^,C'^), given by 

{A»cf>){x):=Acf>{R^'x) . (2) 

Then, on the space of 4-spinors L'^{R^, C*) = L^{R^, C^) © L'^{R^, C^), one can define 
the following unitary representation (denoted again by •) 



(a.(^)](x) - A^«^)(^A_/^^<^(^I'^) 



(3) 



The radial symmetry of V implies that H'^ + V commutes with •. The eigenspaces AA^ 
are thus SU{2) invariant. Now, let J = {Ji, J2-, J3) be the total angular momentum 
operator associated to the representation •. The eigenvalues of = Jf + J2 + Ji are 
the numbers (j^ — 1/4) , where j takes all positive integer values. If ip is an eigenvector 
of with eigenvalue (j^ — 1 /4) , then the SU (2) orbit of </> generates an SU (2) invariant 
complex subspace of dimension 2j > 2. This implies the following fact, which will be 
used repeatedly in the present paper: 

Lemma 1 // </> G L'^(R^,C'^) is not the zero function, then there is A £ SU{2) such 
that (j) and A* cj) are two linearly independent functions. 

Proof of the Lemma. Assume, by contradiction, that C(/> is SU{2) invariant. Then 
is an eigenvector of Ji for i — 1,2,3, hence it is eigenvector of J^. But we have 
seen that in such a case, the SU{2) orbit of </> must contain at least two independent 
vectors: this is absurd. □ 

As a consequence of the Lemma, the spaces Ml have complex dimension at least 
2. The degeneracy is higher in general: for each j > 1 , -\- V has infinitely many 
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eigenvalues of multiplicity at least 2j. Note that in the case of the Coulomb potential, 
the eigenvalues are even more degenerate (see e.g. |16p. Now, on the Grassmannian 

manifold 

Gn{H^'^) := {W subspace of H^'^ {W^ ,£%dimc [W) = TV} 
we define the Dirac-Fock energy £^ as follows 

N 



^=l 



(4) 



2 J JR3xE3 |a; - y\ 

where k > is a small constant, equal to e^/47reo in our system of units, {ipi, ■ ■ ■ V'Jv} 
is any orthonormal basis of W, ^ denotes the iV-uple {ipi, . . . V'iv), P* is a scalar and 
Rq, is a 4 X 4 complex matrix, given by 

N N 
P*W = {'^ei^)^'^ii^)) ' R-fi^iV) = Y i'eix) '^■<P*iy) ■ (5) 

i=i e=i 

Saying that the basis {ipi, • ■ • i/'w} is orthonormal is equivalent to saying that 

Grami2^' = l^v ■ (6) 

We will use interchangeably the notations £^{W) or £^{'^). The energy can be 
considered as a function of W only, because if m G U{N) is a unitary matrix, 

£^,{u^) = £',{^) . (7) 

with the notation {u'^)k — Ukiipi- 

Note that since V is radial, the DF functional is also invariant under the represen- 
tation • defined above. Its set of critical points will thus be a union of SU{2) orbits. 

Finally let us introduce a set of projectors as follows: 

Definition 2 Let P he an orthogonal projector in L^(M'^,C'*), whose restriction to 
is a hounded operator on H ^ (R'^, C*). Given e > 0, P is said to he e-close 
to A'j^ := X(o,+oo)(^^'') if and only if, for all ip e H^M.^,C^), 



(~c'A + c^Y(p-Al]^ <s (-c^A + c^Yi; 



L2(R3_C4) 

In IT^ the following result is proved 



L2(R3,C4) 



Theorem 3 (||14|) Take V , N fixed. For c large and eq, k small enough, for all P 
Eq-cIosc to A^, 

c{P) := inf sup £':{W) 

is independent of P and we denote it by E'^. Moreover, E'^ is achieved hy a solution 
Wk. =span{'ijji, . . . iPn} of the Dirac-Fock equations: 

H"n,wJ^ = ' < e,^ < 1, 

Grami2vl/ = If, ^ ' 

with 

Hl^^:=(H^ + V + .p^.l^)^^. [ ^IpyMyldy. (MF) 

m Jr3 \x-y\ 
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Remark. It is easy to verify that eo > given, for c large and k small enough, 
X{a,oo){H^^wJ is Eo-closeto A+. 

Corollary 4 (^14j) Take V, N fixed. Choose c large and k small enough. If we define 
the projector 

with given by formula ifM_F)) . then 

El^ m\n £1{W)^ min St{W) . (2) 

J=+ W = W ^ solution of (DF) 

Another variational problem was introduced in the works of Bach et al. and Bar- 
baroux et al. (HnHnilll) : define 

^« = K.w = ^[^^oo){Hly^) ; W e Gn{H^'^)}, (3) 

and 



Sl~ := {7 e Si{L^) ,7 = 7*, K.w^ ^ ' 

P+~ 7P"~ = , -P^~ < 7 < P+~ , tr 7 = TV}, 

with the notation P~~ := I — , and 5*1 being the Banach space of trace-class 
operators on L'^(M.^, C^). For all 7 G S^~, let 

Here, Pj{x) := 7^,5(2;, 2;) = Y,n ''^n\ipn{x)\'^ , with w„ the eigenvalues of 7 and V'n 

the eigenspinors of 7, and 7(3;, y) = J2n 'Wn4'n{x) eg) ipn{y), i.e., 7(x, y) is the kernel of 
7- 

In [TH it has been proved that for every P'^— G T'k , the infimum of JFf on the set 



More precisely, with the above notations. 



iS* ~ is actually equal to the infimum defined in the smaller class of Slater determinants. 

K..W 



Theorem 5 ( jl2p For k small enough and for all P~^^ G 'Pk, one has 

inf inf £^,iW) (4) 

Moreover, the infimum is achieved by a solution of the projected Dirac-Fock equations, 
namely 

N 
i=l 

with P^^i^i = V'i (i^ I,. ■ ■ N), and for Wmin := span{il>i, . . . , ■f/'Ar) , 



Grami2\I' = In 



(5) 



Let us now define the following sup-inf: 



sup inf ei{W). (6) 



Then, Theorem |31 has the following consequence: 
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Corollary 6 If k is small enough, 



el = sup inf T^{-f) 

From the above definitions, Theorem |3| CoroUary 0] and the remark made after 
Theorem O we clearly see that for all k small and c large, 

K > <■ (7) 

One can hope more: 

Conjecture M: The energy levels and coincide, and there is a solution of 
the DF equations such that 

£l{W^.)^e%= inf £1{V). 

VeGN(P+„,eHi/2) 

In other words, the max-min level ejj is attained by a pair {W, P^^) such that W = W. 

This paper is devoted to discussing this conjecture, which, if it were true, would 
allow us to interpret the Dirac-Fock model as a variational approximation of QED. In 

order to study the different cases that can appear when studying the problems and 
ej: for k small, we begin by discussing the case k = 0. 

Proposition 7 Conjecture M is true in the case k — 0. 

Proof. The case k = is obvious. Indeed, all projectors coincide with the 

projector X[o.oo){H'^ + V). The level Eg, seen as the minimum of Corollary 2, is 
achieved by any A^-dimensional space Wmin spanned by TV orthogonal eigenvectors 
of H'^ + V whose eigenvalues are the N first positive eigenvalues of H'^ + V , counted 
with multiplicity. Then E^ is the sum of these N first positive eigenvalues. Clearly, 
{Wrmn,X[o,oo){H'' + V)) realizes eg. □ 

The interesting case is, of course, k > , when electronic interaction is taken into 
account. For At > and small two very different situations occur, depending on the 
number TV of electrons. 

The first situation (perturbation from the linear closed shell atom) corresponds 

to 

I 

N = ^N^, I eZ+ (8) 

is treated in detail in §3. 

We recall that iVf is the dimension of the eigenspace Ai'i — Ker(i/'^ + — 1) 
already defined. Under assumption ©, for k = 0, there is a unique solution, W§, to 
the variational problems defining Eq and eg, 

/ 

i=l 

The "shells" of energy , !<«</, are "closed": each one is occupied by the 
maximal number of electrons allowed by the Pauli exclusion principle. The subspace 
Wq is invariant under the representation • of SU{2). 

We are interested in solutions of the Dirac-Fock equations lying in a neighbor- 
hood C Gn{H^^'^) of W§, for K small. Using the implicit function theorem, we are 
going to show that for each k small, W^: exists, is unique, and is a smooth function of 

K. 

Information about the properties enjoyed by is given by 
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Proposition 8 Fix c large enough. Under assumption for k small enough, 

inf S'.iW), (9) 

weGjv(P+„,-effi/2) 

and is the unique solution of this minimization problem. 

This proposition will be proved in §3. Our first main result follows from it : 

Theorem 9 Under assumption for c > fixed and k small enough, = and 
both variational problems are achieved by the same solution of the self- consistent 
Dirac-Fock equations. For e^, the optimal projector in Vk, is P^w"- 

Proof. The above proposition implies that for k small, 

e«> inf e^iW) = £^,{W^J ^ E^. (10) 

Therefore, ej: — E^. Moreover, by Proposition |H1 is achieved by a couple {W^,P) 
such that P = P^-^^c , being a solution of the Dirac-Fock equations. This ends the 
proof. □ 

The second situation (perturbation from the linear open shell case) occurs when 

TV = ^ iVf + A:, / e Z+, < fc < N^i^^ . (11) 
1=1 

It is treated in detail in §4. 

When (|llll holds and when k = 0, there exists a manifold of solutions, 5*0, whose 
elements are the spaces 

1=1 

for all Wfj^^ G Gt{Mjj^-y). These spaces are all the solutions of the variational prob- 
lems defining Eq and eg. The (/-I- l)-th "shell" of energy X'jj^i is "open": it is occupied 
by k electrons, while the Pauli exclusion principle would allow Nf^-^ — k more. JVote 
that we use the expression "open shell" in the linear case k = only : indeed, adapt- 
ing an idea of Bach et al. |17| . one can easily see that for k positive and small, the 
solutions to (DF) at the minimal level E^ have no unfilled shells. 

For K > and small we look for solutions of the DF equations near Sq (see §4) . We 
could simply quote the existence results of and show the convergence of solutions 
of (DF) at level E^, towards points of Sq, as k goes to 0. But we prefer to give another 
existence proof, using tools from bifurcation theory. This approach gives a more precise 
picture of the set of solutions to (DF) near the level E^ f Theorem 112(1 . 

In particular, we obtain in this way all the solutions of (DF) with smallest energy 
E^ (Proposition inj). 

We now choose one of these minimizers, and we call it W^. We have -P~^c (W^) = . 
Since V is radial, belongs to an SU{2) orbit of minimizers. We are interested in 
cases where this orbit is not reduced to a point. Then the mean-field operator 
should not commute with the action • of SU (2), and one expects the following property 
to hold: 

(P) .■ Given c large enough, if k is small, then for any solution of (DF) at level 
E^, there is a matrix A G SU{2) such that 

P^.yvM'W^J^O. (12) 
Let us explain why (P) contradicts Conjecture M: 
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Proposition 10 // (P) is satisfied, then for c large enough and k small, given any 
solution of the nonlinear Dirac-Fock equations such that £^{W^) — E^, we have 

E^^^£^,{W^)> inf f^(W^). (13) 



This proposition will be proved in §4. Moreover, we verify (see Proposition ll5() that 
(P) holds when / > 1 and k — 1, i.e. when in the linear case there is a single electron 
in the highest nonempty shell. 

Our second main result follows directly from Propositions II (Jl and 1151 

Theorem 11 Take 

I 

1=1 

For c large and k > small, there is no solution of the nonlinear Dirac-Fock 
equations with positive Lagrange multipliers, such that the couple 

realizes the max-min e^. 5*0 Conjecture M is wrong. 



3 - Perturbation from the Unear closed shells case. 

Let us recall that we are in the case 

iV = ^7Vf, /ez+, 

being the dimension of the eigenspace = Ker(iJ'^ + l^ — 1). We want to apply 
the implicit function theorem in a neighborhood of W§ , for k small. For this purpose, 
we need a local chart near W§. Take an orthonormal basis (V'l, • • • , iPn) of Wq, whose 
elements are eigenvectors of H'^ + V, the associated eigenvalues being < ■ ■ ■ < fiN 
(i.e. A^, . . . , AJ counted with multiplicity). Let Z be the orthogonal space of Wq for 
the L2 scalar product, in H^^^{R^,C^). Then Z is a Hilbert space for the H^^^ scalar 
product. The map 

C : X = (Xi, ■ ■ • ,Xn) span{iji +Xi, ' ■ ' ,i^N + Xn) , 

defined on a small neighborhood O of in Z^, is the desired local chart. Denote 
the N X N matrix of scalar products {xiiXi)l'^- Then 

The differential of this functional defines a smooth map F^, : O d Z^ [Z')^ , 
where Z' C H^^/"^ is the topological dual of Z for the H^/^ topology, identified with 
the orthogonal space of for the duality product in H~^/^ x H^^'^. Note that F„ 
depends smoothly on the parameter k. A subspace C(x) is solution of (DF) if and only 
if Ff^(x) = 0. To apply the implicit function theorem, we just have to check that the 
operator L := Dy^Fo{0) is an isomorphism from to its dual (Z')^. This operator 
is simply the Hessian of the DF energy expressed in our local coordinates: 

Lx = {{H, + V- ^ll)xu ■■■ AHc + V- ^in)xn) ■ (14) 

Under assumption the scalars fik, k = 1, . . . N, are not eigenvalues of the restriction 
of H'^ + V to the L^-orthogonal subspace of W§. This implies that L is an isomorphism. 
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As a consequence, there exists a neighborhood of W§ x {0} in Gn{H^^^) x R, x 
(— Ko,Ko) and a smooth function h'^ : {—ko,ko) ~^ ^ such that for k G (— ko,ko), 
:= /i'^(k) is the unique solution of the Dirac-Fock equations in il. Moreover, for all 
K e (— Ko,Ko)i the following holds: 

uiW^) = , Wu e SU{2) . (15) 

Indeed, the subset A of parameters k such that 1)15(1 holds is obviously nonempty (it 
contains 0) and closed in (— ko,ko). Now, for k in a small neighborhood of A, the 
SU{2) orbit of stays in fl. But this orbit consists of solutions of the Dirac-Fock 
equations, so, by uniqueness in fl, it is reduced to a point. This shows that A is also 
open. A is thus the whole interval of parameters (— kq, kq). 

Now we are in the position to prove Proposition |S1 

Proof of Proposition ISl 

Remember that for k = 0, PQ^rc coincides with X(o,oo)(-f^'^ + ^)- Now, Wq is clearly 
the unique minimizer of £g on the Grassmannian submanifold Gq Gn{P(Iw^-^^^^)- 
More precisely, in topological terms, for any neighborhood V of Wq in Gn{H^^^), there 
is a constant 6 = S{V) > such that 

Soiw) > s^iw^) + s ,yw eG+ n {Gn{h^'^) \ v) . (i6) 

Moreover, looking at formula H14|l . one easily sees that the Hessian of £q on G^ is 
positive definite at W§ . We now take k > small, and we consider again the chart 
C constructed above. We define the submanifold := GjqiP^y/aH'^^'^)- Then the 
restriction C'^ of G to [P'^^^Z)'^ is a local chart of near W^. For k. small enough, 
there is a neighborhood U oi Q in such that the second derivative of £^ o G+ is 
positive definite on W+ := U C\ (P^y^^Z)^ . The functional £^ o G+ is thus strictly 
convex on Now, for k small, there is a unique Xk G i^t such that G+(xk) = W^. 
Then the derivative of o G+ vanishes at Xk- As a consequence = G+(xk) is the 
unique minimizer of £^ on :— C'^{U^). Now, we choose, as neighborhood of W§ 
in Gn{H^^^), the set V := G{U), and we consider the constant S > such that itTCjl is 
satisfied. Taking k > even smaller, we can impose 

mm£^^ + d/2< inf 

V+ G+\V+ 

Hence, is the unique solution to the minimization problem □ 

4 - Bifurcation from the linear open shell case. 

Recall that here we are in the case 

TV = ^ iV^ + A:, / e Z+, < fc < iVp+i . 
For K — 0, there exists a manifold of solutions, Sq, whose elements are the spaces 

for all Wj^-^ 1^ G Gi{Mf_^_i). These spaces are all the solutions of the variational prob- 
lems defining Eq and Cq. 

For K > and small we want to find solutions of the DF equations near 5*0, by using 
tools from bifurcation theory. 
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If A/4.1 has only multiplicity 2, then implies k = 1 and by Lemma 1 of §2, Sq 
is an SU{2) orbit. Then, as in §3, one can find, in a neighborhood of Sq, a unique 
SU{2) orbit 5^ of solutions of (DF). But there are also more degenerate cases in which 
A/4-1 has a higher multiplicity, and Sq contains a continuum of SU(2) orbits. In such 
situations, k = is a bifurcation point, and one expects, according to bifurcation 
theory, that the manifold of solutions Sq will break up for k ^ 0, and that there will 
only remain a finite number of SU(2) orbits of solutions. To find these orbits, one 
usually starts with a Lyapunov-Schmidt reduction: one builds a suitable manifold 
which is diffeomorphic to Sq (see e.g. JH|)- When Sq contains several SU{2) orbits, 
the points of S'k are not necessarily solutions of (DF), but contains all the solutions 
sufficiently close to 5o. Moreover, all critical points of the restriction of £^ to 5"^ are 
solutions of (DF). The submanifold 5^ is constructed thanks to the implicit function 
theorem. More precisely, we consider the projector 11 : ^ ®{=i -^i ■ To each 
point z e S'o we associate the submanifold := {w € Gn{H^^^) : Hw — z}. For w a 
point of F^, let A„, T^F^ C T^Gn{H'^''^). Then the following holds: 

Theorem 12 Under the above assumptions, there exist a neighborhood of Sq in 
G]\[{H^^^), a small constant kq > 0, and a smooth function h : Sq x (—kq,kq) 
such that 

(a) h{z,0) = z Vz e S'o 

(b) Denoting S'k := h{So, k), S'^ is also the set of all points w in U, such that 

<(£:^)'(H,^>=0, VCeA^ (17) 

(c) h{z,K) e F^ , V(z,k) e Sq X {-kq,kq). 
Proof. 

We first fix a point z in Sq. Let M be the orthogonal space of 0[^i in H^/"^ for 
the scalar product. As in §3, we can define a local chart Cz ■ O C {JV)^ Fz near 
z, by the formula C{x) = spanijp + x), where ip = ("01; ' ' ' j V'w) is an orthonormal basis 
of z consisting of eigenvectors of H'^ + V, with eigenvalues ^1 < • ■ ■ < hn (j'-e. A^, . . . , Af 
counted with multiplicity). The Hessian of £§ o Cz at x = is given once again by 
formula (|14(l . It is an isomorphism between (JV)^ and its dual. So, arguing as in §3, 
we find, by the implicit function theorem, a small constant > 0, a neighborhood lUz 
of z in Fz and a function hz : {—Kz,Kz)—^^z such that: 

(i) hM^z 

(ii) hz{K) is the unique point w in fJ^ such that 

< (£^)'(ii;),e>=0, VeeA^ (18) 

Since 5*0 is compact and £^{w) a smooth function of (w,k), it is possible to choose 
Kz , ^z such that kq := inizeSa '^z > 0, with f7 := Y^jzi^Sa ^ neighborhood of Sq, and 
h{z,K) :— hz{n) a smooth function on Sq x {—kq,kq) with values in Q. This function 
satisfies (a,b,c). □ 

From (b) any critical point of £^ in f2 must lie on S^. From (c) it follows that S'„ 
is a submanifold diffeomorphic to Sq, and transverse to each fiber Fz in Giq{H^/^). 
If z e 5*0 is a critical point of o h[-, k), then, taking w = h{z, k), the derivative of 
£^ at w vanishes on T^S^.. From (b), it also vanishes on the subspace A^, which is 
transverse to TyjS^, in TzGn{H^^'^)^ hence {£'^)'{w) = 0. This shows that the set of 
critical points of £^ in coincides with the set of critical points of the restriction of £^ 
to S'k. Arguing as in the proof of Proposition [SJ one gets more: 

Proposition 13 For k > small, the solutions of (DF) of smallest energy are 
exactly the minimizers of £^ on S^, • 
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We are now ready to prove PropositionllUI 

Proof of Proposition 1101 . Since k is small, for any matrix A e SU (2) the map 

P/^AmW induces a diffeoniorphism between the submanifolds GNiP^W"^^^^) ^^'^ 

Now, we fix A G SU{2) such that (|12|l holds. Then there exists a unique point 
W+ G Gjv(i?^/^) such that 

PZw^W+ = 0, P+A,wsW+ ^A.W^ (19) 

By (IT^ . we have 

W+ ^A» Wl- 

On the other hand, in HI it was proved that 

Sl{A,Wl)= sup Zl{W) (20) 

and A • is the unique solution of this maximization problem. Therefore, 

El{A.W'i,)>£l{W^). 

But 

£l[W^) > inf £^,{W) , 

weGN(P+„,,ffi/2) 

hence, by invariance of £^ under the action of SU{2), 

E^, = £^,{A.W^)> inf £^,{W), 

and the Proposition is proved. □ 

Since there are no solutions of (DF) under level E^, and ej: < E^, PropositionllUI 
has the following consequence: 

Corollary 14 // (P) is satisfied, then for c large enough and n small, there is no 
solution W, of the nonlinear Dirac-Fock equations with positive Lagrange multipliers, 
such that the couple 

iw*,p;t.wj 

realizes the max-min e^. So Conjecture M is wrong when (P) holds. 

We now exhibit a case where (P) holds. 

Proposition 15 Assume that N = ^N^ + 1, I >1. Then (P) is satisfied. 

1=1 

Proof. 

Step 0. Fix c large enough and take a sequence of positive parameters (k£)£>o con- 
verging to 0. Let {W^)e>o be a sequence in GN{Et'^^^), with a minimizer of £^^ on 
. Let -0^ e W[ be an eigenvector of the mean-field Hamiltonian H^^ , normalised 
in and corresponding to the highest occupied level. Extracting a subsequence if nec- 
essary, we may assume that ip^ ^ ip'^ £ ■^'i+i — Ker(i?° + V — Aj_,_j). Moreover, from 
Theorem II 21 we have 



^ Wo" = A^? © C . 
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Step 1. Fix c > 1 . Since P^^ vF'^V'l = 0) we can write, by a classical result due to 
Kato, 



1 



+ 00 



(21) 



1 f^°° 

= ^ f°^iH' + V- iri)-\^A,ws - ^ws){H' + V - ir?)-^" dr^ + o(k^) , 
where by Q.w we denote the nonlinear part of H^ -^: 

Hl^ = H'' + V + Knw. 
But note that since the space 0{=i Alf is invariant under the action of SU{2), 

So, wc just have to prove that for c sufficiently large and for all tp"^ € Mj_^_-^ , there 
exists A G SU{2) such that 

/ + 00 
{H' + V- i7i)'\nA.^^ - ^^^){H' + V- irj)-^ij' drj^O. (22) 
-OO 



Since 



{H^ + V - ijj)-'^'' = 



•^7+1 - iV 



and fl^o = > 



what we need to prove is that for all nonzero ip'^ € Mj^-^ , there exists A € SU{2) such 
that £°(^1a»v= ''P^) 7^ 0) with 



/ + 00 
-OO 



-1 dri 



•^7+1 - in 



Step 2. We give an asymptotic expression for when c — > +oo: 



-^/_:(>'-)-^)-^4( 



where L^, in the Fourier domain, is the operator of multiplication by the matrix 

/ + 00 
{-iu + 13+ {a- p)/c)-\-iu + 1)"^ du . 
-OO 

Here, we have used the standard fact that 

A' 



(23) 



(24) 



We have 



(-m + /3+(a-p)/c)-^ 



—iu + W^(p) —iu — u)'^{p) 



A"-(P) 



with 



.^^(p):=v/rTwW, Ac^(,) = --(^)^(^+^--^)/--) . 
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Hence, by the residues theorem, 



TT 



step 3. It is well known (see ^Hl) that ■0'^ can be written as 



2c 

4> € L^(R'^, C^) being an eigenstate of {^ + V), with eigenvalue = linic^+oo {^"i+i — 
c^). Since we have assumed that V is smooth, this asymptotic result holds for the 
topology of the Schwartz space 5(R^). So, 

where 

fiA, 0) {\A . 0p * ^) - (< A . 0, >c. (A . 0) . (25) 

What remains to prove is : 

Step 4. For any eigenvector 4> of the Schrodinger operator —^ + V, there exists an 
A £ SU{2) such that f{A, (f>) ^ . 

Proof of Step 4. We consider the integral 

lA,<t,{r) ■■= < {x- a)(j), f{A, (j)) >c2 {ruj)duj . 



Since (f> has exponential fall-off at infinity, the electrostatic field jA^^p * takes the 

asymptotic form ^ J^^ 1^" '^1^) ]^ "I" '^(jip') when \x\ is large. The same phenomenon 
holds for the convolution product < Au (pjcj) >c2 ■ As a consequence, for r large. 



^ < ^ • 0, (/) >C2 ) ( y ^ < 0, ^ • >C2 {ruj) du?j 



Since • is unitary, the Cauchy-Schwartz inequality gives 
By Lemma 1 of §1, we can choose A such that <j) and A» (p are not colinear. Then 



\4)\ {r Lo) duj ^ / \A u 4>\ {r Lu) duj > j < A • (j), 4> {r uj) duj 
S2 Js^ 



So there is a constant 5 > such that, for r large enough, 

\rlAAr)\>s(f \^?)(! \<t>?{roj)d,j). (26) 

Being an eigenvector of the Schrodinger operator —^ + V^ the function cannot 
have compact support. So the lower estimate H26I) implies that the function IaA''') 
not identically 0, hence f{A, (j>) ^ . Step 4 is thus proved, and (P) is satisfied. □ 

Aknowledgement. The authors wish to thank the referee for useful comments on 
the first version of this paper. 
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